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AN  EXACT  LOWER  CONFIDENCE  BOUND  FOR  THE  RELIABILITY  OF  A  SERIES  SYSTEM 
WHERE  EACH  COMPONENT  HAS  AN  EXPONENTIAL  TIME  TO 
FAILURE  DISTRIBUTION 

by 

Tapas  K.  Sarkar 

1.  Introduction  and  Summary 

Let  us  consider  a  complex  system  consisting  of  k  components.  The 

components  are  connected  in  series  so  that  the  system  fails  if  at  least 

th 

one  component  fails.  Let  T^  represent  the  time  to  failure  of  the  1 
component.  The  reliability  of  the  system  at  time  t  is  given  by 

(1)  R(t)  =  P[TX  >  t,T2  >  t,...,Tk  >  t]. 

Biis  Implies  that  T,  the  time  to  failure  of  the  system,  is  given  by 

(2)  T  =  min  [ T .  ]  and  R(t)  =  P[T>t]. 

i=l,...,k 

let  us  assume  that  T^  has  an  exponential  distribution  with  parameter 
so  that 

t 

(3)  P(Ti  >  t)  =  e  ,  i  =  1,.  ..,k. 

Initially,  let  us  assume  that  are  independently  distributed. 

Then  from  (l)  and  (2), 

-(\  +. .  .+X,  )t 

(b)  R(t)  =  P[T  >  t]  =  e  1  , 


so  that  the  time  to  failure  of  the  system  has  an  exponential  distribution 
with  parameter  ^  +  '  *  '  +  In  general,  will  not  be 


known.  We  want  to  find  a  lower  confidence  bound  for  R(t). 

We  assume  that  the  system  as  a  whole  is  either  impossible  or  too 
expensive  to  test.  Hence,  no  data  is  available  for  T.  However,  one  can 
test  each  individual  component  and  observations  are  available  for  T.  for 
each  i  =  l,...,k.  We  want  our  confidence  bound  to  be  based  on  these  latter 
observations. 

Kraemer  [7]  has  given  two  methods  for  obtaining  the  lower  confidence 
bound  of  the  above  reliability.  However,  her  methods  do  not  give  an  exact 
lower  confidence  bound.  The  method  proposed  by  us  will  give  an  exact  lower 
confidence  bound.  It  will  also  be  shown  that  our  method  works  under  cer¬ 
tain  circumstances  even  when  are  not  independently  distributed. 

We  have  given  numerical  examples  to  compare  Kraemer' s  first  method 
with  ours.  We  have  also  compared  our  lower  confidence  bound  with  that 
given  by  maximum-likelihood  estimates. 

Kraemer' s  second  method  is  based  on  a  version  of  Tchebyshev's  in¬ 
equality.  It  appears  to  us  that  in  most  cases  it  is  inferior  to  Kraemer 's 
first  method.  Accordingly,  we  will  not  mention  it  any  further. 

Lentner  and  Buehler  [8]  have  used  the  Lehmann-Scheffe  theory  to  ob¬ 
tain  a  uniformly  most  powerful  unbiased  test  for  the  hypothesis 
0.  Inverting  the  test  they  obtained  an  exact  confidence  bound  for  the 
reliability  of  the  system.  The  method  has  been  extended  to  the  case  of 
k  components  by  El  Mawaziny  [1].  Unfortunately,  the  methods  are  com¬ 
plicated  and  even  in  the  case  of  k  -  2  and  very  small  sample  sizes  one 
has  to  solve  equations  by  trial  and  error  to  obtain  the  confidence  bounds 
for  0.  For  larger  values  of  k  and/or  for  moderate  sample  sizes 
elaborate  computer  programming  is  required.  In  fact,  El  Mawaziny  and 
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Buehler  [2]  derived  an  approximate  method  when  there  are  k(k  >  2) 
components  because  of  the  difficulty  in  computation  for  their  exact  meth' 
ods.  It  will  be  shown  below  that  our  method  is  very  easy  to  compute  and 
can  be  easily  worked  out  with  a  desk  calculator.  The  methods  of  [8]  and 
[1]  are  applicable  only  when  the  components  are  independent.  As  will  be 
shown,  our  method  also  works  for  dependent  components. 

El  Mawaziny  and  Buehler  have  used  a  theorem  of  R.  von  Mises  to 
obtain  large  sample  approximations  of  the  methods  of  [8]  and  [1].  Their 
method  depends  on  asymptotic  expansions  of  certain  quantities.  It  turns 
out  that  if  the  first  two  terms  of  the  expansions  are  used  (and  the 
higher  order  terms  are  ignored),  the  confidence  bound  of  reliability  is 
the  same  as  that  given  by  the  maximum-likelihood  method.  The  approxima¬ 
tion  can  be  improved  by  using  higher  order  terms,  but  the  degree  of 
improvement  is  not  known. 

2.  The  Method 

To  illustrate  the  method,  let  us  assume  that  k  =  2. 

Let  T^  be  independent  exponentially  distributed  random  variables 
with  parameters  i  =  1,2.  Let  T  =  min(T^,Tg).  Then 

-(\,+\2  )t 

P(T  >  t)  =  P(TX  >  t,  T2  >  t)  =  e 

Hence,  T  is  exponentially  distributed  with  parameter  (k^+X^). 

Let  us  assume  that  n^  similar  components  of  type  i  are  on  test 
simultaneously.  The  test  is  without  replacement  so  that  as  one  compo¬ 
nent  fails  it  is  removed  from  test  and  is  not  replaced.  We  terminate 
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the  test  as  soon  as  among  n^  have  failed.  This  method  is  known  in 
the  literature  as  type  II  censoring  of  r^  out  of  n^  without  replace- 


Let  X/,  v  <  <  . . .  <  X/  x  be  the  order  statistics  of  the  times 

(1)  (2)  (i^) 

to  failure  of  the  components  of  the  first  type,  and  Y^  <  Y^j  <  ••• 

<  Y/  x  be  the  same  for  the  components  of  the  second  type.  For  our 
v  *2 ' 

method  we  require  that  ri  =  r2‘  Let  r  =  min(r^,rg).  Hence,  we  work 
with  termination  points  r  for  each  component.  Kraemer's  method  or 
the  methods  of  [1],  [2],  and  [8],  however,  do  not  require  this  restric¬ 
tion.  For  simplicity,  let  us  assume  that  nj_  =  n2  =  n.  As  will  be  evi¬ 
dent,  there  is  no  loss  of  generality  in  making  this  second  assumption. 


Di  -  xar 

dj  =  x(j)  '  Vi)’  J  = 


Bj  =  (n-J+l)Dj  ,  j  =  l,...,r. 


Epstein  and  Sobel  [3],  and  Epstein  [4]  have  shown  that  B^,...,B  are 
independently,  identically  distributed.  The  distribution  is  exponential 


with  parameter  Similarly,  let 


E1  =  Y(l)’ 


Ej  =  Y(J)  '  VD'  J  =2 . r> 

=  (n-j+l)Ej  ,  j  =  l,...,r. 


Tlien  , . . . , Tr  are  independently  exponentially  distributed  with  parameter 


i 

i 


mm##*? 


Tl,”‘,Tr  are  ■'•n^ePen<^en't 


\g*  Also, 

Let 

(7)  Zj  =  min(6j,Tj). 


h’hen,  as  we  have  already  stated,  z  ,  .,.,zr  are  independently  exponen¬ 
tially  distributed  with  parameter  Let 

r 

( 8 )  rV  =  £  z  . . 

J=1  J 

It  is  well  known  that  2rV(\^+\2)  has  chi-square  distribution  with  2r 

2 

degrees  of  freedom.  Let  x  denote  a  random  variable  distributed  as  chi- 

2 

square  with  2r  degrees  of  freedom.  Let  y^^a)  be  such  that 

(9)  Ptx2  <  x|r(a)]  =  l-a. 


Hence, 


or 


P[2rV(\1+\2)  <  )|r(a)]  =  l-a 


x|r(«) 

P[A^+\2  ^  2rV  ^  =  ^-Q! 


-(\l+\2)t 

and,  since  R(t)  =  e  is  a  decreasing  function  of  +  a. 2  for 

fixed  t,  we  get 


(10) 


-(X^+Xgjt 
P[e  >  exp 


l-a. 


This  gives  an  exact  lOO(l-a)  per  cent  lower  confidence  bound  on  the  reli¬ 
ability  of  the  system  at  time  t  when  the  .system  consists  of  two  compo¬ 
nents  in  series  and  the  time  to  failure  distribution  of  each  component  is 
exponential  and  the  two  exponentials  are  independent. 

It  is  obvious  that  the  above  method  does  not  require  the  assumption 
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are  unequal,  then  in  (5)  and.  (6)  the 


that  n^  = 
definition  of 


(11) 


=  n.  If  and  n^ 

6.  and  t,  will  be  changed  as  follows: 
J  J 


B,  =  (n1-j+l)D. 


T  .  = 


J  \ 


<VJ+1)Ej  1 


j  =  1,2, • • .,r. 


■flie  rest  of  the  theory  goes  through  and  (7),  (8),  (9)>  and  (10)  hold 

with  the  new  definitions  of  & . ,  t.,  and  z.. 

J  J  J 

3.  Extensions  of  the  Method 

A.  When  the  number  of  components  is  larger  than  2. 

Let  us  consider  a  system  consisting  of  k  components  in  series.  Let 
be  the  time  to  failure  of  the  i  component.  Assume  that  T^,Tg,...,T^ 
are  independent  exponentially  distributed  random  variables  having  parameters 
\1,\2,...,\k  respectively. 

Let  T  =  min(T^ , . .  .  ,Tk ) .  Then 

P(T  >  t)  n  P(TX>  t,...,Tk  >  t). 

Since  j  •  •  •  ,T^  are  independent, 

P(T  >  t)  =  P(T1  >  t)  ...  P(Tk  >  t) 

-(\^+  ...  +  A^Jt 

-  e 


Hence,  T  is  exponentially  distributed  with  parameter  +  ...  +  s  \, 
say. 

Suppose  we  have  n  components  of  each  type  on  test,  and  let  us 
perform  type  II  censoring  of  r  out  of  n  without  replacement.  Let 
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"Hy,"  wnmvmmmtLmmn.m..,, 


our  data  consist  of  <  xj^j  < 

tion  2,  for  i  =  l,2,...,k,  let 


.  <  X/^l,  i  =  As  in  sec¬ 

tor 


D(i)  _  X(D 

ui  ~  h 


(12) 


D(i)  -  X<iJ  -  X(i)  1-2  r 

J  ‘  (J)  (j-l)’  J  1 


(13) 


6^  =  (n-j+ljD^  ,  j  =  l,...,r, 

z  =  rain  6^  ,  j  =  1, . . .  ,r. 
J  i=l,...,k  J 


z,  ,...,z  are  independently  exponentially  distributed  with  parameter  \. 
l  r 


Let 

(14) 


rV  =  Sz  . 
i=l 


Bien  2rV\  has  chi-square  distribution  with  2r  degrees  of  freedom,  and 


(15) 


P[e  exp  \  -t 


•  -1 
i 


.4^1 

2rV  j 


]  =  1-a. 


Biis  gives  us  an  exact  lOO(l-a)  per  cent  lower  confidence  bound  for  the 
reliability  of  this  system  at  time  t. 

If  n^,...,n^  are  different,  then  the  modification  proposed  in  the 
last  paragraph  of  section  2  applies. 


B.  When  k  =  2  and  the  distributions  are  not  independent. 

Marshall  and  Olkin  [9]  and  Harris  [6]  have  considered  the  following 
model  of  a  bivariate  exponential  distribution. 

Let  us  consider  3  independent  series  of  events  happening  in  time. 

Let  us  consider  the  series  i,  i  =  1,2,3.  Let  represent  the  interval 

between  the  n^*1  and  the  (n-l)^  event  of  type  i.  We  assume  that  llfn^ 
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has  exponential  distribution  with  parameter  X^,  independent  of  n.  With 

an  arbitrary  time  origin  let  us  assume  that  the  time  to  the  first 

event  of  type  i,  has  the  same  distribution  as  mentioned  above.  Further, 

/  \ 

let  us  assume  that  for  fixed  i,  the  distributions  of  '  are  independ¬ 
ent  for  ail  n.  Hence,  we  have  three  independent  Poisson  series  of  events 

and  U.  has  distribution  exp(X. ) ,  i  =  1,2,3  and  U. ,  IL,  U  are 
l  1  i  c.  2 

mutually  independent. 

Let  us  now  consider  a  system  with  two  components  in  series.  Let  us 
assume  that  component  1  fails  if  either  or  happens,  and  compo¬ 
nent  2  fails  if  either  Ug  or  happens.  The  system  fails  if  at  least 

one  component  fails. 


Let  T^,  i  =  1,2,  be  the  time  to  failure  of  the  i  component. 
Then  ^  =  min(U-L,U^),  T2  =  min(U2,U^)  and 

(16)  P[TX  >  s,  T2  >  t]  =  P[UX>  s,  U2>  t,  U3>max(s,t)] 

=  exp[-X1s  -  Xgt  -  X^  max(s,t)]. 


This  defines  the  distribution  of  a  bivariate  exponential  distribution. 
The  reliability  of  the  system  at  time  t  is 


(17)  R(t)  =  P(TX>  t,  T2>  t)  =  exp[-(X1  +  X2  +  X  )t]. 

This  shows  that  the  time  to  failure  of  the  system  has  exponential  distribu¬ 
tion  with  parameter  X^  +  X2  +  X^.  Also,  the  marginal  distributions  of 
T^  and  T2  are  exponential  with  parameters  X^  +  X^  and  X2  +  X^, 
respectively.  T^  and  T2  are  not  independently  distributed. 

Let  T  =  min(T^,T2).  Then 

(18)  P(T  >  t)  =  P(TX>  t,  T2>  t)  =  exp[-(X1  +  X2  +  X  )t], 
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(18)  is  really  the  same  as  (17). 

Now  suppose  we  have  the  same  data  as  in  section  2.  We  define  , 
t  ,  Zy  and.  rV  as  in  (5),  (6),  (7)>  and  (9).  Then  Zy...,zr  are 
independently  exponentially  distributed  with  parameter  X^  +  X^  +  X^. 

With  this  modification  we  get  that  2rV(\1  +  X2  +  X^)  has  chi-square 
distribution  with  2r  degrees  of  freedom  and,  as  in  (9)  and  (10),  we  get 


(19) 


P[e 


-(X^Xg+X^Jt 


>  exp 


-t 


Xgr<a> 

2rV 


|  ]  =  1-a. 


ftiis  gives  an  exact  lOO(l-a)  per  cent  lower  confidence  bound  for  R(t). 
When  n^  and  ^  are  different,  we  proceed  as  in  the  last  part  of 
section  2. 


C.  As  in  case  B,  but  k  >  2. 

Hie  bivariate  exponential  distribution  has  been  generalized  by 
Marshall  and  Olkin  [9]  to  multivariate  exponential  distribution. 

Let  us  consider  a  system  with  k  components  in  series.  Let  S 
denote  the  set  of  vectors  (s^,...,s^),  where  each  Sj  =  0  or  1,  but 
(s^,...^)  f  (0,...,0).  For  any  vector  s  e  S,  max(x^  3^)  is  the 
maximum  of  the  x^'s  for  which  s^  =  1.  The  multivariate  exponential 
distribution  is  defined  as 

(20)  P(T1  >  x^...,^  >  x^)  =  exp[-  E  Xs  max(xi  s^]. 

Here  T^,...,T^  denote  the  time  to  failure  of  components  l,...,k.  The 
marginal  distribution  of  any  subset  of  , . . . , T^  is  again  a  multivariate 
exponential  distribution  and,  in  particular,  tne  one-dimensional  marginal 
distributions  are  exponential  and  are  not  independent.  As  in  case  B, 
we  can  generate  the  distribution  in  (20)  by  considering  2  -1  independent 
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series  of  Poisson  events  and  the  associated  failures  of  components.  For 
example,  for  three  components  we  need  7  series.  Hie  associated  failures 
are 


1 

2 

3 

4 

5 

6 

7 


shock 

1 

2 

3 

4 

5 

6 

7 


component 

1 

2 

3 


1,2 


1.3 

2.3 


1,2,3 


parame ter 


The  reliability  of  the  system  at  time  t  is  given  by 


(21) 

As  in  case  A, 
j  =  1, • • • ,r, 


R( t)  =  exp[ -  E  \  .  t). 
seS  S 

if  we  have  data  for  the  k  components  and  define 

i  =  l,...,k,  then  z  =  min  j  =  1,.. 

J  i=l, . . . ,k  J 


,r,  has 


exponential  distribution  with  parameter  \  =  E 

seS 

proposed  under  case  A,  section  3,  Is  applicable. 


\ 

s 


and  the  method 


D.  When  type  II  censoring  with  replacement. 

Let  us  suppose  that  n  components  of  a  type  are  on  test  and  we 

replace  each  failed  component  immediately  by  a  new  one.  We  stop  as  soon 

as  r  failures  are  obtained.  This  is  termed  in  the  literature  as  type  II 

censoring  of  r  out  of  n  with  replacement.  Let  D,  be  defined  as  in 

J 

(5)-  Let  us  define 

(22)  =  nDj,  J  =  l,...,r. 

10 


(•WfF- 


t  f*rr  'i 


Epstein  [4]  has  shown  that  6^,...,6  are  independently  exponentially 

distributed  with  the  same  distribution  as  the  parent  population.  With 

this  new  definition  of  6.  (and  t.,  the  methods  of  sections  2  and 

J  J  J 

3  go  through  and  we  can  obtain  an  exact  lOO(l-a)  per  cent  lower  confidence 
bound  for  the  reliability  of  the  system.  As  before,  the  method  works  when 

nlJ,’*,nk  are  n°^  etlua^• 

4.  Confidence  Bounds  Based  on  Maximum  Likelihood  Estimates 

We  will  consider  case  A  of  section  3  and  for  illustration  we  will 

take  k  =  2. 

Let  us  consider  component  i  and  let  the  order  statistics  obtained 
by  type  II  censoring  of  r  out  of  n  without  replacement  be  denoted  by 

v(^ )  <  v(l)  <  <  v(i)  ,  _  ,  O 

X(i)<x(2)<  ...  <X(r),  1  -1,2. 

Using  the  transformation  given  in  (5),  the  likelihood  of  the  sample  is 
proportional  to 

L  =  exp[-xi(&[l)  +  ...  +  6^)],  i  =  1,2, 

where  6^  ,...,6^^  are  independently  identically  distributed  as 
exponential  with  parameter  k^.  Hence,  the  maximum  likelihood  estimate 
of  k^  is  given  by 


6  =  J  6^)  ,  x  =  E  5^ ' , 
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Then  6  and  t  have  gamma  distributions  with  density  functions 

..  r  Rr-1  -Xt6 

X  o  e  1 

f(&)  = - (r-i)i  "  >  &  >  0,  X2  >  0, 


(25) 


Now, 


\r  Tr*1e"^T 

g(T)  -  '"(r-l): - ’  1  i.0’  X2>0' 


3  log  L  _  _r_  .  _  ,  p 

2  -  -  ,  l  -  i,d. 

a  k  k 


—  is,  in  a  sense,  the  "asymptotic  variance"  of  X. .  However,  since 
1  o  *  d  1 


xt  \ 

X^  is  not  known,  we  replace  —  by  — 


Definition:  Let  be  a  random  variable  depending  on  n  and  let 
X  be  any  other  random  variable.  Hien  we  define 


X  as  n  -»  » 


if,  as  n  -»  co,  the  cumulative  distribution  function  of  X  tends  to  the 

n 

cumulative  distribution  function  of  X  at  the  points  of  continuity  of 

P 

the  latter  function.  Also,  we  define  X^  -+  X  if,  given  any  e  >  0, 


P(  |x  -x|  >  e)  -»  0  as  n  -»  ». 


We  shall  prove 


Theorem  i.l: 

Let  6^,...,B^  and  -^,...,1  be  two  independent  random  samples 
from  exponential  distributions  with  parameters  X^  and  respectively. 

Then 
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1m,i>  v: immmmmtrnt***™  ** +%*-*** wy t~ <*"•**■  r  t 


|  +  \  -  \  -  X 

T  =- - = - —  ^  N(0,l)  as  r 

r  77 - ~ 

J  2  +  2 

V  6  T 


r  (1)  r  (2) 
where  6  =  E  '  and  t  =  E  & .  ' . 

i=l  i=l 

We  will  require  the  following  theorem  to  prove  Theorem  4.1: 

Theorem  4.2: 

Suppose  5  =  l,...,r)  is  a  sample  from  a  k-dimensional 

distribution  with  finite  means  j  and  positive  definite  covariance 
matrix  H^jH,  i,j  =  l,...,k.  Let  g(x1,...,xk)  be  a  function  which 
possesses  first  derivatives  =  g^ ,  say,  i  =  l,...,k,  at  all  points 

in  some  neighborhood  of  (^,...,4^),  and  let  g.^  =  g^  .  . . ,  . 

Then  if  at  least  one  of  the  g°  is  /  0, 

'/r[g(x],,...,xk)  -  gC^,-..,^)] 


X  .,rn  „  oo, 

-»  N[0,  E  a,  .g.  g ,  J  as  r  -»  ®. 

i,j=l  J  J 


Proof:  See  [l.'j,P-  260,  section  9- 3 -la]  • 


Proof  of  Theorem  4.1: 

Let  &  =  &/r,  t  =  r/r.  From  what  we  are  given, 


E<V  -  >  v<5i)  -  4  > 


E(u)  -  V(  t  )  =  4-  . 


y  "  K 


J  "  x.2 
X2 


Cov(&.,t  )  =  Oj  i,j  =  l,...,r. 
^  J 
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Consider  g(x,y) 


Then 


and 


o 

S1 


2J/2  ^  °’  - 


/.2  .2,1/2  ^  °> 


+  ®. 


Hence,  by  Theorem  4.2,  as  r  -»  ®, 


Now 


v/r[g(6,t)  -  g(i  ,  ^  )] 


^  N[0,  2  cr  g°g°]. 

i,  j=l  J  J 


g(5^)  = 


r  r  .  . 

£  +  7  '  VX2 


/  2  2 
/r  r 

I*  2  +  2 

5  T 


and 


g(r  >  r  )  =  °> 


V  ^2 


2 

£  o. 


o  o 


1,J=1  1J  1  J  ^ 


Hence, 


g( 5 ^ T )  -  g(r  >b  )] 
\  X2 


i 


r  r 

-2  +  2 
O  T 


^  N( 0,1 ) , 


as  r  -+  ®, 


4 


' "»  » .«wf  trwwf.wTPfoswiwM.* 


and.  the  theorem  is  proved. 

Theorem  4.1  gives  us  an  asymptotically  exact  confidence  interval  for  tne 
reliability  of  the  system.  Let  k^  be  such  that 

j'k“  JL. 

-  ^ 


Then 


-  k  ]  =  1-a 
cr 


) 


(26)  P[\1+X2  <  \  +  \  +  k aJ^2  +  -2]  =  l-o, 

(27)  P[e'<klVt>  exp{-t(§  *  f  +  k,  ^TT^,)}  1  -  l-«- 

(27)  gives  the  lOO(l-o)  per  cent  lower  confidence  bound  for  the  reliability. 
Hie  bound  is  exact  if  r  -»  °°.  The  method  of  this  section  generalizes  when 
the  number  of  components  in  series  is  more  than  2. 

Since  this  confidence  interval  (  h.wnr  ionfidence  bound,  l)  in 
based  on  the  maximum  likelihood  estimate,  for  large  samples  it  is  the 
smallest  confidence  interval.  It  is  not  known  how  good  this  bound  is  for 
small  samples. 

We  will  make  come  comparison  of  this  method  with  the  method  of 
section  2. 

When  the  number  of  items  on  test  for  each  component  is  different, 
the  remark  at  the  end  of  section  2  applies. 


5.  Kraemer's  Method  [7] 

For  the  case  A  of  section  Kraemer's  first  method  is  as  follows: 


15 


2r.  a. 
i  1 


is  distributed  as  chi-square  with 


2r. 


1 


degrees  of  freedom. 


Let  A  =  min  (r.a.).  Then 

1=1,2,..., k 


-(\1+...+Lk)t 

P[e  >  exp 


-t 


2 

X2(r1+...«k)(a) 

2A 


|]  >  1-Of, 


where,  if  U  is  a  chi-square  variable  with  2(r^+  .. 

2 

of  freedom,  then  Xo/  \(or)  defined  as 

e-\ r2+  •  •  •  +  1  ^  / 


+  r^)  degrees 


P[U  < 


+  rk} 


(<*)]  =  1-cr. 


6.  A  Comparison  of  the  Methods 

We  will  compare  the  methods  of  sections  2,  4,  and  5  for  a  system  with 
two  independent  components  in  series.  The  time  to  failure  of  each  compo¬ 
nent  is  distributed  exponentially  with  parameter  i  =  1,2.  Ihe 

comparison  will  be  based  on  the  expected  values  of  the  upper  confidence 
bounds  for  under  each  method.  Since  the  left-hand  end  point  is 

zero,  this  gives  the  expected  length  of  the  confidence  interval  for  + 
Since  the  reliability  is  a  (strictly)  monotone  function  of 
this  furnishes  a  valid  comparison. 


Exact  Method  of  Section  2 


x|r(o) 

P[\l  +  \2  <  2ry  ] 


1-a. 
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4C¥*  WT f  . 


Since  z^  are  independent,  identical  exponential  variables  with  param¬ 


eter 


^1  +  ^2*  (8).  rV  has  gamma  distribution  and 


(28) 


where 


^rX2r(cy)i  Xgr(»)  !  i  X2r(a)  Xi  +  X2 
2rV  ^  "  2  ^  =  2 


I  z4 


2  z, 

i=i 


r  -  1 


„  ,  \r  r-1  '(WX 

„/  1  v  i«  1  ^1+X2^  x 

E(  — 5  =  Jo  *  • - n^TTi - to 


r-l  r-2  (^1+^2  )x 


^1  +  ^2  p®  (X-^Xg)  x 

r-l  Jo 


(r-2)! 


dx 


h  +  x2 

r-l 


,  provided  r  >  2. 


Kraemer's  Method  of  Section  5 


where 


X4r(<*) 

+  \2  <  i  >  l-«> 


A  =  min(ra1,  ra2) 


■  min[2^;  •  2r“iV  2 ^  •  2r°2x2J- 


2 


Since  2ra^\^  are  distributed  as  \  with  2r  degrees  of  freedom  and 
are  independent. for  i  =  1,2, 


P(A  >  x)  =  P(2ra1\1  >  2^x)  P(2ra2\2  >  2\2x) 


■  'J^x  ur'le'UduHJT2x  v  e 


Hence,  the  probability  density  function  of  A  is  given  by 


dv  ] . 
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I 

I 


,  -X.,  x  ,  -\_V 

,r  r-1  1  ,r  r-1  2 

x  e  e 

f(x) - (r-l):  Jx  (r-1)!  dv 


,  -X.0X  V 

. r  r-1  2  .  r  r-1  1 

X2X  e  «  xiv  e 

*  — ( r-i ) ;  Jx  (r-1):  dv'  x  i  °- 


Integrating  by  parts, 


,  x 

, r  r-1  1  .  .  /x  \S 

X.,x  e  r-1  -k  x  (\0x) 

f<x>  -  (r-1):  z  e  -7T- 

S=0 


1  -\«x 

,  r  r-1  2  .  .  /.  v£ 

e  r-1  -X^x  (X^x) 

+  — TTTF-  s  e 

s=o 


s; 


Hence, 


(29)  E(i)  =  £  |  f(x)dx  =  2  for  r  >  2, 

A  Jo  x  s-o  ( r-i ) !  s.  ^  f+s-1 


and 


. <r<“\  <r<«>  (r«-2>:  Kte  *  X2X1 

(30)  Et-gr- )  -  -5—  jrti)!':  ^  r«-l’  proVlded 


r  >  2 


s=o 


Asymptotic  Method  of  Section  4 


nr-v  ,  ,  r  r  .  r  r  •> 

P[X1  +  X2  <  5  +  7  +  ka  J~2+  ~2] 


=  l-a> 


Using  (25), 

(31) 


E(l  +  7J  =  (xi +  M' 


Substituting  5  =  P  cos  9,  1  =  p  sin  9,  and  using  (25)1  we  get 

Vb' 


E[  /-2+  -2]  =  E[ 
8  t 


y/r  Vb2  +  T2 


5T 


Vr  (^2)  (2r-2)!  w/2  (cos  Q  sin  Q) 

=  (r-i):  (r-i):  J.  : 


r-2 


d9 


\2r-l 


(\^  cos  9  +  \  sin  9)£ 

provided  r  >  2. 
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I 

I 

I 

a 

a 


(Put  t  =  tan  £•) 


y?  (\1\2)r(2r-2): 
"Tr-1)!  (r-1)! 


2r-1  pi  tr~2(i-t2)r~2(ut2)2dt 
)o  (a+t)2r'1-(b-t)2r-1 


,2r 


- 1; 


where 


k0  O  k— 

a  =  /l  +  (A  -  4  >  °, 
X1  *1 


I 

I 

a 

a 

B 

Q 

a 

o 

o 

Q 

0 

n 


.  .  ,X 2v2  X2  .  . 

bN1  +  (  h 


We  will  break  the  integrand  into  partial  fractions: 

t^rtTW)2  _  2r_1 


or 


2r-l  B. 


=  £ 


+  £ 


(a+t)2r-1(b-t)2r"1  i=l  (a+t)1  j»l  (b-t)J 


rz  (r'2)  (-l)sttr4£s'2  +  2tr+2s  +  tr+2s+2J 
s 

s=o 


=  £  A  (a+t)2r"1'i(b-t)2r’1 

t=l 


2r-l 

1  _  /v  .»2r-l-J/  ,  v2r-l 

+  £  B.(b-t)  (a+t)  , 

J-l  J 


where  A' a  and  B's  are  independent  of  t. 

We  have  Ur-2  unknowns.  We  will  differentiate  n  times 
(n  x  0,1, . . . ,2r-2 )  with  respect  to  t  and  put  t  =  -a  and  t  =  b. 

This  will  give  us  two  systems  of  2r-l  equations  each.  Each  system  is 
triangular  and  we  can  easily  solve  sequentially  for  Agr  ^,Bg r  ^,. . . ,A^,B^. 
For  the  n  differentiation  the  general  pair  of  equations  is 


f 

t 
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r-2 

£ 

s=o 


z  (r;2)  (-D3  [(r+f-2)(-a)»2s-2-" 


and 


♦  2(r+2s)(-a)r+2s-n  ♦  (rt2s+2)(-a)rt2s+2-n] 


k=o 


r* -P 

„  ^r-2^  ^ js^r+2s-2j^r+2s-2-n 

s=o  3 


~,r+2s v,  r+2s-n  , r+2s+2 v,  r+2s+2-n , 

+  2<  n  >b  +  <  n  )h  ) 


"  /2r-lw  N2r-l-(n-k) /  vk 

=„ E  B2r-l-k  <  n-k><atb>  (-1) 

k=o 


with  the  proviso  that  (^)  =0  if  p  <  q,.  After  solving  for  A's 
we  get,  using  ( 31 ) , 


<*=>  *$*$**'&$! 


-jr  (X.  \ )r(2r-2)!  a1"1 
*  +  “a  - * 


(r-l)l  (r-l)! 


2r-l 


t2^'1  ITT  <4^  -  -7^ *  A!  <¥> 

1=2  a  (a+1) 

2r-J.  Bj  /  1  _  1  \  ,  b  N 1 

+  T,  1-1  i-1  1-1^"  ®l  b-1  ^ 

+  L  J_1  (b-i)J  1  bJ  1  1  b  1 


j=2 


for  r  >  2. 


We  will  now  make  some  numerical  comparisons  between  (28),  (30) 
(32). 
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and  B's, 


,  and 


nr  i 


Numerical  Comparisons 

Uie  following  table  gives  some  numerical  values  of  E  (upper  99  Per 
cent  confidence  bound  for  X^  +  Xg)  based  on  (28),  (30),  and  (32)-  In 
(30)  we  have  used  tables  of  cumulative  binomial  distribution  a:  follows: 


B(n,k,p)  =  £  (g)pS(l-p)n"S- 

s=k  s 

r'1  (r+s-2).'  X1  x2  +  X2X1 

Z  ( r-I ) fs'i  ~  ■  Tr+s-'l’"  '  r  ^  2 

s=o  v  '  (Xx+X2) 


r-1  0  X,  ^  X0 

1  v  /r+s-2 w _ l_\r-l  /  2_vS 

"  r-1  r-2  AX1+Xp;  '  ( X1 +X^ 


X_  r-1  _  X_  ,  Xn 

-  *  (”s22)  <rf\r>  (rjr)s 

s=o  12  12 


k  <2r2><^>5  <^>2r-2-8 


*  k  2\  <2r;2>^>8  <^>2r-2-s 


X  X 

=  — 7  B(2r-2,r-l,  - — 7-) 
r-1  v  ’  ’  X^+X^ 

+  ^  li.BU/.a,/,  A-)]  . 

In  (32),  instead  of  using  the  partial  fraction  expansion  we  have  performed 
numerical  integration  of 


j>n/2  (cos  9  sin  6)r  2d© 

°  (X1  cos  ©  +  Xg  sin  ©)2r  ^ 


We  have  used  Simpson's  rule  after  subdividing  the  interval  [0,ji/2] 
into  ten  equal  parts. 


\2  r 

E  (exact) 

E  (Kraemer) 

E  (asymptoti 

2 

2.5  5 

10.30 

11-37 

8.72 

10 

7.85 

8.40 

6.90 

20 

6.6  0 

6.98 

* 

30 

6.1b 

6.48 

* 

1 

5  5 

13-73 

19-66 

12.20 

10 

10.47 

15-49 

10.47 

20 

8.80 

13.41 

* 

30 

8.18 

12.64 

* 

4 

4  5 

4.584 

5.004 

3.874 

10 

3.494 

3-674 

3.064 

20 

2.944 

3.034 

# 

30 

2-734 

2.794 

* 

* 

For  larger  values  of  r  the  error  involved  in  approximate  integra¬ 
tion  becomes  large  and  hence  these  values  have  not  been  given. 

As  has  been  pointed  out  before,  the  above  table  also  gives  a  valid 
comparison  between  the  confidence  intervals  of  R(t)  as  given  by  the 
different  methods. 

In  the  range  of  a  that  we  are  interested  in  (e.g.,  1-ce  =  .90, 

•95,  etc),  it  is  expected  that  the  exact  method  will  give  a  shorter 
confidence  interval  than  the  Kraemer  method.  Since  Kraemer's  method  is  an 
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i 

0 

4 

approximate  one,  in  general,  for  X.  and  \  widely  different  (as  in 
the  second  set  of  examples),  the  exact  method  is  expected  to  be  signif¬ 
icantly  superior  to  Kraeraer's.  This  happens  because  the  error  of  ap¬ 
proximation  in  her  method  exceeds  the  precision  obtained  by  using  a 
larger  number  of  degrees  of  freedom  in  the  chi-square  distribution. 

When  and  are  close  (as  in  the  first  set  of  examples),  some¬ 

times  Kraemer's  method  may  give  a  better  confidence  interval  than  the 
exact  method;  but  the  difference  between  the  two  methods  is  not  expected 
to  be  large.  For  very  small  values  of  a  (e.g.,  1-a  =  .99,  -999>  etc  ), 
there  is  some  numerical  evidence  that  sometimes  Kraemer's  method  is 
slightly  better  than  the  exact  method  when  \  -  It  may  be  pointed 

out  that  for  =  \g,the  difference  between  the  expected  value  of  the 
upper  confidence  bound  for  given  by  Kraemer's  method,  and  the 

expected  value  as  given  by  the  exact  method,  goes  to  zero  as  r  goes  to 
infinity. 

Kraemer's  method  has  one  advantage  over  the  exact  method.  The  for¬ 
mer  is  applicable  even  when  the  termination  points  (n)  for  the  compo¬ 
nents  are  different.  For  the  exact  method  some  data  have  to  be  ignored 
to  equalize  the  termination  points. 

For  large  r  it  is  expected  that  the  asymptotic  method  will  be  bet¬ 
ter  than  the  exact  method.  Without  a  knowledge  of  the  rate  of  convergence 
of  the  exact  and  asymptotic  methods,  it  is  not  possible  to  assert  to  any 
degree  of  precision  how  the  two  methods  will  compare.  However,  it  is 
expected  that  the  exact  method  will  compare  favorably  with  the  asymptotic 
method  for  large  r. 

Since  the  methods  of  [1]  and  [8]  are  based  on  inverting  uniformly 
most  powerful  unbiased  tests,  they  may  give  better  lower  confidence  bound.- 

£ 

4 
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than  those  given  by  the  method  proposed  in  this  paper.  This  is  especially 
true  for  small  values  of  r.  For  large  values  of  r  the  confidence  bounds 
given  in  [1]  and  [8]  will  tend  to  the  bounds  given  by  the  method  of  [2], 
and  will  be  the  same  as  those  given  by  the  asymptotic  method.  As  we  have 
already  pointed  out,  the  methods  of  [1]  and  [8]  are  computationally 
complicated  even  for  the  case  of  two  components,  and  work  only  when  the 
components  are  independent.  In  fact,  the  computational  difficulties 
motivated  the  authors  to  get  (approximately)  asymptotic  results  in  their 
subsequent  joint  paper  [2].  These  three  methods,  however,  do  not  require 
the  same  termination  points  for  every  component. 

When  V,  -»  0,  the  exact  method  and  the  asymptotic  method  proposed 
in  this  paper  can  be  easily  compared.  If  A. 2  -»  0,  then  we  are  reduced 
to  the  case  of  a  single  component  with  exponential  time  to  failure  distribu¬ 
tion.  In  this  case,  ®  for  all  i,  and 


(10)  reduces  to 


z. 

l 

rV 


min(6.,T.)  =  &i 
r 

E  8.  . 

,  ,  i 


Xpr(«) 

P[e  >  exp(-t-  2-—  )]  =  1-oc  , 


the  well-known  result  about  a  single  component. 

Let  0  i-n  (27)-  Then 

P[e  1  >  exp  £-t(|  +  k^^|)  ^  ]  =  l-o. 

Following  Halperin  [5], 

v/?(R(t)  -  &(t))ou£  N(0,1). 


2h 
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Hence,  as  r  -»  », 

exp j-t 
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proposed  by  Lentner,  Buebler,  and  El  Mawaziny. 
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